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Abstract. 

Charged polymer brushes are layers of surface-tethered chains. Experimental systems are 
frequently strongly charged. Here we calculate phase diagrams for such brushes in terms of salt 
concentration n s , grafting density and polymer backbone charge density. Electrostatic stiffen- 
ing and counterion condensation effects arise which are absent from weakly charged brushes. 
In various phases chains are locally or globally fully stretched and brush height H has unique 
scaling forms; at higher salt concentrations we find H ~ n s , in good agreement with exper- 
iment. 



Introduction. - Charged polymer layers are a major concern of polymer science, featuring 
in numerous applications such as colloid stabilization, surface modification technologies, mem- 
brane preparation [1] and emerging biotechnologies [2] . A much-studied class of layer is the 
polymer brush, an assembly of chains end-tethered to a substrate (see fig- - Charged brushes 
are important in biosensor technologies such as DNA microarrays [2] and have provided model 
systems to study collective physics of charged polymers at interfaces [3-9] . Properties such as 
brush height have been measured by the surface force apparatus, neutron scattering and other 
techniques [3-5] . The principal theoretical frameworks have been self consistent field theory 
and scaling theory [6]. Scaling theories identified distinct brush "phases" [7,8] and later a full 
"phase" diagram of brush types [9] . 

Most of this theoretical work has addressed weakly charged brushes where local Gaussian 
chain statistics are only weakly perturbed by electrical forces. Real polymers, however, are 
frequently strongly charged. The natural measure of this is the reduced backbone density 
("Manning parameter") qo = where I is the (monovalent) charge spacing and the Bjerrum 
length Zb = e 2 /e w kT = 7 A in water. Strongly charged polymers include single-stranded 
DNA (q w 1.6 [10]) and the widely studied polystyrensulfonate, PSS (e.g. q w 2.8 at 80% 
sulfonation [4]). When go is not small, two qualitatively new effects arise, (i) Strong repulsive 
forces may stretch chains beyond the linear Gaussian regime into rodlike configurations with 
end-to-end size of order the chain contour length itself: the polymer size saturates, (ii) 
Manning condensation [11]. When qo > 1, a rodlikc polymer attracts its own counterions 
so strongly that a fraction condenses into a region close to the polymer, renormalizing the 
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effective Manning parameter to unity, qo — ► 1. The implications of (i) and (ii) have been 
explored theoretically in certain brush regimes [12-14]. 

The aim of this letter is to establish complete phase diagrams for strongly charged brushes. 
One of our motivations is the urgent need to establish physics underlying DNA microarrays, 
a technology based on DNA layers of enormous importance in health care, drug development 
and basic biological research. We will show that the saturation and Manning condensation 
effects described above, (i) and (ii), lead to brush phases whose scaling forms for brush height 
H and counterion layer thickness D have no counterparts among weakly charged brushes [9] . 
Two brush classes are identified, defined by the magnitude of the intrinsic persistence length, 
a. Polymers in the "strong" class (a > ^b) ar e strongly perturbed by charges, suffering length 
saturation for high enough backbone charge density whereas for "weak" systems (a < l-g) 
Manning condensation intervenes before saturation can onset. Correspondingly, each class has 
a distinct brush phase diagram. The key point is that typical experimental brush systems are 
"strong" and their chains locally fully stretched, invalidating results based on flexible Gaussian 
polymer concepts. Brush height then depends on the salt-dependent persistence length l p for 
a semi-flexible polymer. The dependence of / p on salt concentration n s remains a highly 
controversial question [15,16]. Here we take l p w £, a relation with considerable empirical 
support [10,17], where £ = (87tZb«s) is the Debye screening length. For the particularly 
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experimentally important "quasineutral" brushes our theory then predicts H ~ n s , in 
good agreement with experiment [3,4]. Interestingly, no such agreement is obtained if one 
uses the often cited alternative quadratic OSF form [15], l p ~ £ 2 . Thus brush experimental 
data comes down rather firmly on the side of the l p ~ £ relation. 

Single Chains. - To appreciate saturation and condensation effects, it is helpful to first 
review results for a single charged chain of contour length L and charge Q — L/l [18]. Through- 
out, we assume ideal Gaussian chain statistics in the absence of charges. For small charge, 
the chain size is Rq « (La) 1 / 2 . When Q becomes large, internal repulsions exert a stretching 
force /i n t ~ IbQ 2 /H 2 [9] (setting kT — 1) where H is the extended chain size. Balancing this 
with the elastic restoring force f e \ « H/Rq and using /mt^b ~ 1 gives the electrostatic "blob" 
size 

r h » a(ZA sat ) 2/3 , U = (M 1/2 • (1) 

For small scales r < statistics are essentially Gaussian, while on larger scales the polymer 
is electrostatically stretched into a linear string of blobs of length [18] 

Hnb*(l*/l) 2 / 3 ; b*=La/l B , l* = l%/a. (2) 

However, this result for H cannot be correct for all charge densities: from eq. Q when 
I = isat the blob shrinks to one persistence length a and the chain size becomes L. This is size 
saturation; further charging produces no further extension. (In this study, the window between 
onset of non-Gaussian behavior and complete size saturation is ignored for simplicity.) 

This suggests H rs L for I > Z sa t< However, this reasoning neglects Manning condensation. 
Since the charge per blob is ~ r\lal, the string of blobs has an effective linear charge density 
coarse-grained over the blob size equal to 

qf « (r/0 1/3 , r = 1%/a . (3) 

Counterion condensation is triggered when = 1, i.e., when 1 = 1*. If we now consider a 
chain charging process (decreasing I from I = oo), this leads us to define 2 classes of polymer: 
(1) The weak class, a < Iq. In this case I* > l sat , so Manning condensation onsets before 
saturation. Thus eq. (J2J for chain size remains valid until I = I* when the polymer reaches its 
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Fig. 1 - Brush phase diagrams for no added salt (grafting density b , backbone charge separation 
I, logarithmic axes), (a) Weak and (b) strong systems. OB: osmotic brush; PB: Pincus brush; 
IC: independent chains; MC: Manning condensation; PMC: partial MC; FR: fully relaxed chains. 
Saturated phases are indicated by prefix "Sat". 



maximum size, b* . Charging beyond I* generates just enough condensate to pin the effective 
charge spacing to I*. H is pinned to b* and saturation never occurs. (2) The strong class, 
a > Iq. This includes most polymer species, e.g. PSS (a ps 12 A) [17] and B-DNA (a ps 450 
A) [10]. Saturation now occurs before condensation can intervene, since I* < i S at- Thus eq. 
J2J is valid until saturation at I = l sat . For higher charge densities, H ps L. 

We remark that the treatment of Manning condensation for weak systems above is oversim- 
plified. One can go beyond this by treating the string of blobs as a cylinder of radius tv Adapt- 
ing Ramanathan's expression [19] for the potential at a charged cylinder, ip(qQ S ) = 2q\n(^/r\y) 
with q — min(l, <?g ff ), a charging process then yields the free energy F = (2Q — H/Iq) ln(£/rt>) 
for qg S > 1. Taking the derivative gives the stretching force /i n t and balancing this with f e \ 
we find H = b* ln(£/rt,), confirming that apart from the log factor chain size is indeed pinned 
to 6*. Such logarithms, which arise elsewhere {e.g. eq. are ignored in the following for 
simplicity. 

Brushes: Weak Systems. Consider now a brush of chains of the weak class, end- 
tethered to a surface at density b~ 2 . We begin with no added salt. Our results are shown 
as a phase diagram in fig. [Ja). The upper portion, Z > Z*, is the weakly charged region 
previously established by Pincus and Borisov et al. [7-9]. We now review their results. At low 
grafting densities lies the Independent Chains (IC) phase where layer height H\q is essentially 
the single chain result, eq. ©, and the height D of the counterion layer ensuring overall 
neutrality is the Gouy-Chapman length Aqc = b 2 /2ttIbQ- These counterions exert a weak 
external stretching force / ox t ps IbQ 2 /b 2 per polymer, less than the internal force /i„t- To the 
left lies the Pincus brush (PB) where instead / oxt > /int- The IC/PB boundary is / ex t ~ /int, 
or 1 1 1* pa (b* /fe) 3 / 2 . In the Osmotic brush (OB), the high grafting density localizes counterions 
in the brush [9], unlike PB and IC where D > H. Thus the PB/OB boundary is Aqc ~ #pb 
or l/l* PS (6*/6) 4 / 3 [9]. All 3 phases coincide at the triple point P* = (b*,l*). What happens 
in the strongly charged region below I*? The weakly charged brush theory naturally raises a 
question it cannot answer. 

The single chain analysis suggests the answer: the dependence on / is removed by Manning 
condensation. We find this is essentially correct, using similar free energy calculations to 
those described for single chains: for I < I*, brush height is fixed to b* (to within log factors) 
independent of I or b. This is the Manning condensation (MC) phase in fig. ^( a )- It nas 2 
sub-regions where D ps H (for b < b*) and D ps \ gc (for b > b*). At the lowest grafting 
densities the MC phase gives way to the partial Manning condensation (PMC) region; here 
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Fig. 2 - Brush phase diagram, strong systems with salt. The I < Ib region (not shown) maps onto 
I — Ib- Arrows indicate how boundaries evolve with increasing salt concentration, (a) L > £ > Ro, 
(b) £ < Ro. Salt-dominated regions (shaded) are blown up at right. QB (Quasineutral brush) and 
QIC (quasineutral independent) are subdivided into bands a, (3, 7. Since < -Ro, 7 does not exist 
for £ > Ro (see text). 



the counterions are so weakly attracted to the surface (D ^> H) that condensate evaporation 
onsets [20]. We find onset at b w 6*/e 1 / 2 , where the small parameter e = 1 / ln(b* / Ib) ■ Deep 
into the PMC region, enough condensate is lost that length saturation occurs, despite this 
being a weak system. Finally, the boundary of the uppermost "fully relaxed" (FR) region is 
defined by H w R . 

Brushes: Strong Systems. - Now we consider strong systems, a > Ib (see fig. H(b)). 
As expected from the single chain discussion, at high charge densities new saturated phases 
appear (indicated by prefix "Sat") where brush height reaches its maximum possible value, 
H rs L. The IC region crosses over to Sat-IC at I « l sat (see eq. (QJ), while since chains in 
the OB have one charge per blob [9] saturation is at Z « a (recall, blob size equals persistence 
length at saturation) giving the Sat-OB phase. Between these two lies the saturated Pincus 
brush (Sat-PB). Its boundary with Sat-OB is Aqc ~ L or 1/1-q » (L/b) 2 ; and with Sat-IC is 
/hit ~ fcxt or b m L. Saturation effects shift the triple point P* to (L,Ib)- 

The essential feature of strong systems is that Manning condensation onsets after chains 
have saturated and are rod-like, when q — 1, i. e. / = Analogously to the I < I* phases 
for weak systems, below / = 1-q lie the MC and PMC phases (fig. ^ where e = l/ln(L/c?) and 
d is polymer thickness.) 

Salt Effects. - For any real brush screening effects of salt must be considered. For brevity, 
we discuss strong systems only and the typical case £ < L. We find Manning condensation 
then renormalizes any / < Ib to I — > Ib with no evaporation effects. Thus fig. [21 shows only the 
portion of the phase diagram above I — Ib] the region below is mapped vertically onto this 
line. Figs.^a) and (b) show the two possible cases £ > i? an d £ < Ro- 

Let us focus on the region affected by salt, shown shaded in fig. [21 (compare to the no- 
salt case, fig-^l. This expands with increasing salt concentration as indicated and is divided 
into two by the vertical line b = £. To its left lie the salt brush (SB) and saturated salt 
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Fig. 3 - The quasineutral brush (QB) behaves as a neutral brush with effective hard core monomer 
size a e ff depending on polymer charge density and salt concentration. There are 3 regimes, giving 3 
bands in the phase diagram, fig.|2] Band a: length saturation, a e e w £. Band (3: string of electrostatic 
blobs, a e ff ~ £• Band 7: overscreening (£ < rb), a c g « r^/t; 2 . 

brush (Sat-SB) phases where interactions are long-ranged and the stretching force is the 
differential osmotic pressure [9] of ions. The SB, established by weakly charged brush theories 
[7-9], saturates (H » L) at small I, crossing over to the Sat-SB phase unique to strongly 
charged brushes. Equating the salt concentration n s to the brush counterion concentration 
n c w Q/Hb 2 defines the salt/osmotic boundaries. Thus the Sat-SB/Sat-OB boundary is 



Moving to the right of b = £ one enters the quasineutral brush (QB) phases. Since the 
screening length is less than the chain separation &, interactions are now short-ranged so chains 
behave as self and mutually avoiding polymers with a certain effective hard core monomer 
size, a e ff, defined as the separation between any two chain segments where their electrical 
interaction energy is order kT. Thus we adapt the Alexander-de Gennes theory of neutral 
brushes [6] viewing each chain as comprising blobs of size b (see fig. |3J|. Depending on /, 
there are three possible forms for a e g: these are the three bands a, (3, 7 in the phase diagram, 
fig- 13 W1 th corresponding structures shown in fig. (1) In the saturated a band, / < Z sat , 
polymers are locally rod-like as discussed for single chains. Flexibility is recovered on scales 
above £ (recall, we take persistence length l v w £). Thus a e g « £ and there are N c g w L/£ 
effective monomers per polymer. (2) The intermediate band, (3. Here the charge density is 
below saturation and locally the polymer a string of electrostatic blobs of size rb (see eq. Q). 
Thus a c ff ~ £ and N e g ~ (Z sa t//) 2 / 3 L/£. (3) In the overscreened 7 band, screening penetrates 
individual blobs, £ < rb, reducing their unscreened electrostatic energy which by definition 
equals kT. Thus a c g exceeds 77,. From a Flory type calculation we find a c g w rj^/f 2 and 
N eS w La/a 2 cS . 

Borrowing the neutral brush result [6], the brush height is -ffq>B ~ N ga^b~ 2 ^ 3 , giving 



This is one of our principal conclusions. The result for band (3 has been derived previously 
by Borisov et al. [9]. Note the identical scaling for (3 and 7 is coincidental. 

Finally, to the right of QB lies the quasineutral independent chains (QIC) region whose 

boundary is Hqb ~ b. The bands a, /?, 7 remain, and chain size is the standard self-avoiding 
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result R w 7V C ^ a c g. This agrees qualitatively with the structure of a single charged chain 
found in ref. [16] though the predicted exponents differ. 

We remark that if we assume the OSF form for persistence length [15], l p « ^b(£/0 2 ; 
we find a qualitatively similar phase diagram to fig. |3 but now: (i) between the Sat-SB and 



I/Ib « (£A) 2 - 
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QB phases a very narrow saturated "Nematic Brush" region appears where rodlike effective 
monomers nematically order and (ii) new bands appear in the quasi-neutral phases. A crucial 
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difference is that the predicted salt dependence in the QB regime is now Hqb ~ n s ■ 

Experiment. We conclude by comparing theory with two experimental studies of 
charged brushes. (1) Tran et al. [3] used neutron scattering to measure monomer density 
profiles of high density brushes (13 ~ Z> ~ 29 A, 100 ~ L ~ 2000 A) of PSS on silicon at 
sulfonation levels 35% to 62%. Using a representative monomer size [21] 2.1 A, this implies 
charge spacings 3.4 ~ Z ~ 6 A. (2) Balastre et al. [4] studied less dense PSS brushes on mica 
(b w 100 A, 1000 ~ L ~ 1500 A, sulfonation 84%-87%, i.e. I w 2.5 A) with the surface force 
apparatus. 

Now PSS is a strong system (a = 12 A) and both experimental systems are strongly 
charged, I < Zb, thus belonging to the line I = Zb in the strong phase diagram, fig.|21 Hence 
at low salt (£ > b) we predict both systems lie in the Sat-OB phase with fully stretched 
chains. Indeed, the observed low salt brush heights are ~ 0.7L and w 0.9L in the Tran et al. 
and Balastre et al. studies, respectively. This is deep into the non-Gaussian regime [12,13], 
confirming the chains are close to size saturation. 

What happens when salt is added? Our theory predicts the system in effect moves along 
Z = Zb in fig. EI passing from the Sat-OB phase (£ 3> b) to the a band of the QB phase (£ <C b) 
where chains are locally saturated. The predicted brush height H is thus initially independent 
of salt concentration n s , then decaying as Hqb ~ n s (eq. @). Something very close to 
these predictions is seen in both experiments: H was roughly constant at low salt followed 
by a rather sharp transition to an approximate power law H ~ with £ = 0.270 observed 
in ref. [3] and ( = 0.30 to 0.33 in ref. [4]. The transition occurred at &/£ as 3 (ref. [4]) and 
6/^w9 (ref. [3]). Since the Tran et al. systems have very small b, the QB is realized only at 
very high salt; we speculate the somewhat lower exponent £ and higher transition &/£ value 
for this data may originate from non-electrostatic interactions (more important at high salt). 

Throughout, we assumed Gaussian statistics, i.e. O solvents. How would different solvent 
conditions affect the strongly charged phases with saturation? (i) Good solvents are expected 
to produce no qualitative change, since these phases are already locally or globally almost 
fully expanded, (ii) Very poor solvents and ion-ion correlation effects have been predicted 
to induce brush collapse [8,22]. For PSS, for example, water is a poor solvent for the PS 
backbone. However, for high charge densities solubility is recovered [3]; indeed, both PSS 
systems above exhibit highly extended brushes at high charge density [3, 4] as discussed. 

In conclusion, we find strongly charged brushes exhibit unique phases where chains are 
globally or locally almost fully extended into rod-like configurations. Systematic experiments 
[3,4] on strongly charged PSS brushes find brush height H as const, at low salt n S) crossing 
over to H ~ with exponents C as 0.3. It is tempting to interpret this in a weakly charged 
brush picture [7-9] as a transition from the osmotic (H ~ n°) to the salt brush (H ~ n s 1 ^ 3 ). 
We find this is incorrect; the transition is in fact from a saturated osmotic brush (if ~ n°) 

— 1/3 

to a quasineutral brush with local saturation, for which we predict H ~ n s provided the 
persistence length scales linearly with Debye length, Z p ~ £ [10,17]. The observed exponent 
is clearly inconsistent with the OSF [15] form Z p ~ £ 2 which would give a stronger decay, 

tt -1/2 

H ~ n s 

* * * 
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